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Abstract 

We extend Rice Formula for a process that is the sum of a smooth process 
and a pure jump process. We obtain formulas for the mean number of both, 
continuous and discontinuous crossings through a fixed level on a compact time 
interval. An application to the study of the behavior of the tail of the distri- 
bution function of the maximum of the process over a compact time interval 
is considered. Further, we give a generalization, to the non-stationary case, of 
Borovkov-Last's Rice Formula for Piecewise Deterministic Markov Processes. 



1 Introduction 

Counting the number of crossings through a fixed level u by a stochastic process 
and studying the behavior of the maximum of the process on a time interval are 
classical, related, problems in Probability Theory. Nevertheless, few is known, 
in general, about the distribution of the random variable number of crossings. 
Therefore, Rice Formulas, which give expressions for the moments of the number 
of crossings, are very important and useful. 

Roughly speaking. Rice Formulas express the moments of the number of 
crossings in an integral form involving the parameters and local properties of the 
process. It is worth to remark that it reduces a global property as the number of 
crossings to local ones, as the finite-dimensional distribution at a finite number 
of points in the parameter space of the process, and its derivatives. 

Different cases of Rice Formula have been known and used for a long time, 
see [12], the first one, due to Rice in 1944 [Ml [17], involve the Gaussian sta- 
tionary case. The major part of the literature is devoted to the Gaussian case, 
and in particular to the stationary case. Extensions for general Gaussian, non- 
Gaussian processes and fields were stated afterwards, see [3l[T0]. Some recent 
extensions to non-Gaussian processes include the Shot Noise processes [5] and 
the Generalized Hyperbolic Process [T]. 
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The applications of Rice Formulas include telecommunications and signal 
processing [1^1 [H], reliability theory in engineering [TB], oceanography (the 
height of sea waves) [H [13] , physics and astronomy (random mechanics) [TT] , 
(the Shot Noise processes) [S] and (microlensing) [T3] , random systems of poly- 
nomial equations [Hill], etc. 

All these works consider smooth processes of class at least . More recently, 
Borovkov and Last in [6] , see also [7] , studied the continuous crossings through 
a fixed level of a class of processes with jumps named Piecewise Deterministic 
Markov Processes. In spite of including jumps, the structure of these processes 
is quite simple in the sense that randomness only is included at the jump 
instants and magnitudes. 

On the present work, we consider a class of stochastic processes with fi- 
nite intensity of jumps and smooth (stochastic) evolution between these jumps 
(excepting Section 4). More precisely, we consider a process X which can be 
written in the form X = Z + J' where Z is a, process with continuously differ- 
entiable paths and is a pure jump process, independent from Z. That is, Z 
describes the continuous evolution of X and J describes the jumps of X. 

Such a process can cross the level u in a continuous way but it also can cross 
u at one of the jumps. Our interest relies on both, continuous and discontinuous 
crossings through the level u. Up to our knowledge, the discontinuous crossings 
have not been considered in the literature. This work began under the initiative 
of Mario Wschebor, and is largely inspired by his fundamental contribution in 
the field, masterly exposed in [3]. 

We obtain formulas for the mean number of continuous and discontinuous 
crossings through a fixed level m G M by such processes on a compact time inter- 
val, and apply them to the study of the distribution function of the maximum 
of the process. We also include a generalization of Borovkov-Last's Rice-type 
formula to the non-stationary case 

The paper is organized as follows. In section 2 we give some definitions, 
introduce the processes we deal with and present the main result. In Section 
3 we present two examples and the application to the study of the maximum 
(on these examples). Section 4 contains the generalization of Borovkov-Last's 
formula. Finally, Section 5 contains the proofs of the main results. 

2 Preliminaries and Main Result 

Let T > and X = {X{t) : t ^ [0,T]) be a stochastic process defined on the 
interval [0, T]. We assume that X can be written in the form X ^ Z + J' where 
Z and J are independent processes on [0,r] as described below. 

We assume some regularity conditions on the process Z, see p], namely: 

Al the paths of Z are , almost surely. 

A2 We assume that the density Pz{t) (x) is jointly continuous for t E [0, T] 
and a: in a neighborhood of u. Further, assume that for every t, t' g [0, T] 
the joint distribution of {Z{t), Z{t')) has a density Pz{t) z(t')('^' '^^^ich 
is continuous w.r.t. t {t',x,x' fixed) and w.r.t. x at u {t,t',x' fixed). 



2 



A3 for every i e [0,r] there exists a continuous version of the conditional 
expectation 

E{Z{t)\X{t) = x), 
for a; in a neighborhood of u. 
AA the modulus of continuity of Z tends to if (5 0: 

w{Z,5):= sup \Z{t) - Z{s)\ ^ Q. 

0<s<t<T,\t~s\<S 

The jump process J — [Jt : t ^ [Oj^]) is based on a general point process 
(T„,^„)„gN on [0, cxd) X M, and is constructed via a family of Markov kernels 
(P^t)i (""i"^) : ?^ € N as follows, see [S]: set tq = and draw according to 
the distribution ttq, then, conditioned on the resulting value of say xq, draw 
Ti with (conditional) distribution 

P(rie- |eo = a;o) = pW(-)- 

Similarly, conditioned on the values of Co:''"ij say xo,ti, draw with distri- 
bution T^^xotxi.')- Then, conditioned on the preceding values and on = xi 
draw T2 with distribution P; . . (•) and so on. Finally, for t„ < t < r„+i let 
ft — n and 

n 
k=Q 

Hence t„ represents the n-th jump instant and ^„ represents the n-th jump 
magnitude (or increment) of the process J^. 

Equivalently, the kernels (-Pj"^) can be used to obtain the actual value of 
the process at the jump instants, in that case we set J(t„) = 

We can identify the marked point process (t„, J(r„)) with the associated 
random counting measure (RCM for short) X]n=o '^(i-,i../(t„)) on [0, T] x R, where 
S(^t,x) is Dirac Delta measure concentrated at {t,x). There exists a predictable 
random process A on [0, c») x R such that: 

E J fd^i = E J fL{dt,dy) 

for any predictable function /; here L stands for the (random) measure induced 
by A. Usually L is called the compensating measure of the RCM /i. 

Under quite general conditions (for example the absolute continuity of the 
kernels and the finiteness of the intensity of jumps), the compensating measure 
L can be written in terms of ordinary Lebesgue integrals. 

We turn now to the definition of a crossing through a fixed level by a function. 
Consider a cadlag, between jumps, function / : [0, T] — ^ M. We say that / 
has a continuous crossing through the level u at s € (0, T) if / is continuous at s, 
f{s) — u and /'(s) ^ 0; and a discontinuous crossing if {f{s~') — u){f{s) — u) < 0. 
If f'{s) > in the continuous case, or f{s^) < u < /(s) in the discontinuous 
one, we say that / has an up-crossing at s, otherwise the crossing is a down- 
crossing. 
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By the preceding assumptions on the process X , we can apply these defini- 
tions to almost all of its paths. Denote (resp. N!^) the (random) number of 
continuous (resp. discontinuous) crossings through the level u by the process 
X on the interval [0, T] and by the total number of crossings. Analogously 
[/„, U!^, [Dm -D^, ^J^) denote respectively the total number, continuous, 
and discontinuous up-crossings (down-crossings). 

It is clear that Nu = + N^, then 

ENu=EN^+EN^. (1) 

The two terms of the r.h.s. of ([T]) are treated separately and by different 
methods. For the continuous crossings we recall that classical Rice Formula is 
based on Kac Counting Formula for the number of crossings of a function, 
and observe that when the function has jumps (and the value u is not one 
of the lateral limits) Kac Formula counts the number of continuous crossings 
through the level, ignoring the discontinuous ones. On the other hand, for the 
discontinuous crossings we use techniques from point processes theory. 

Our main theorem includes the case of non-Gaussian continuous processes 

Z. 

Theorem 1. Let Z and J he two independent processes on [0,T] such that: 

• Z verifies the conditions Al, A2, A3 and AA, 

• J is a pure jump process, as described above, with finite intensity. 

Then, the mean number of continuous and discontinuous crossings through the 
level u by the process X on the interval [0, T] are given, respectively, by: 

EK=/ E{\X{t)\\X{t)^u)px(t){u)dt 

J[0,T] 

= I E{\Z{t)\ \ X{t) ^ u)pxit){u)dt, 

J[Q,T] 

and 

EiV^=E (( l{{X{t-)~u){X{t-)+y-u) <0}L{dt,dy), 

J J[0,T] xR 

where lA is the indicator function of the set A and L is the compensating 
measure of the random counting measure generated on [0,T] x by the jump 
process J . Similar formulas hold for the number of up and down crossings. 

Some remarks are in order. First, note that when the jump process 
vanishes, that is, if J{t) — almost surely for all t e [0, T], Theorem [T] reduces 
to Classical Rice Formula for the process Z. 

Next, observe that the random variable J{t) does not need to have a density 
for each t, but, in case it does we have a more explicit result. 

Corollary 1. If J(t) has a continuous density pjf^i-^{x) fort e [0,T] and x G M, 
we can also write 

EN^^ [ dt [e (\Z{t)\ I Zit) = v) pzit)iv)pjit){u' v)dv. 

J[Q,T] Jr ^ ' 
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Finally, a careful analysis of the proof of Theorem [T] in Section 5 shows 
that the result in Theorem [T] holds true whenever the law of the process Z, 
restricted to the subintervals [r,;,Ti+i], conditioned to the paths of the jump 
process verifies the hypothesis A\ - AA. Besides, A3 can be weakened assuming 
that the product of the conditional expectation and the density oi X{t) (or Z{t) 
in Corollary [ij is continuous. 

We end this section specializing these results to the case where Z is a Gaus- 
sian process, here, the hypothesis of continuity of the densities and of the con- 
ditional expectation may be released since they follow from the conditions of 
non-degeneracy of the distribution of Z{t), t e [0,T], and on the regularity of 
the paths. Besides, the ingredients in the formulas are computable explicitly. 

Corollary 2. Let Z be a Gaussian process with paths such that for every 
t the distribution of Z{t) is non-degenerated, assume further that J' is a pure 
jump process independent from Z with finite intensity of jumps. Then the result 
of Theorem]^ hold true. 

Corollary 3. Under the hypothesis of Corollary\^ if Z has constant variance, 
in particular if it is stationary, then the formula for the continuous crossings 
reduces to: 



Jo 

Here we used the well known fact that for a centered Gaussian process, hav- 
ing constant variance implies the independence of the process and its derivative 
at each point. 

Remark 1. Let us define the density 



Thus, the mean number of continuous crossings through the level u by X is 
proportional to a probability density function evaluated at u. This expression is 
similar to Rice 's original result, but it involves a different density function. 

Furthermore, observe that p is a mixture of the density of X{t) fort £ [0,T]. 
In particular, if X is a stationary process, the density p reduces to that of X{0), 
just as in the original formula due to Rice. 



3 Examples 

In this section we present two examples of the computation of the formulas in 
Theorem [T] and their application to the study of the distribution of the maxi- 
mum of these processes in a compact time interval. Furthermore, we compare 
which kind of crossings predominate as the level u tends to infinity. 





Then, we can write 



TE\Z{0)\p{u). 
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3.1 Stationeiry processes with 1-dimensional Gaussian dis- 
tribution 

In this example we assume that Z = {Z{t) : t G [0,T]) is a stationary, centered 
Gaussian process with paths. Let T be the covariance function of Z, i.e: 
r(T) = EZ(0)Z(t), assume also that r(0) = 1/2. 

Let i7 = {J{t) : t € [0, T]) be a pure jump process constructed in tlie 
following way. Given a sequence of independent and identically distributed 
random variables : n € N) with common centered normal distribution with 
variance 1/2, and p gM. such that \p\ < 1, define 

^0 = Co, and for n > 1 : A„ = pAn-i + p^ ^n- 

It is easy to see that the sequence (A„ : n G N) is centered, stationary Gaussian 
with variance 1/2 and covariance between and A„ given by p"'/2. 

Further, consider a Poisson process v = {ut •.t> Q) with intensity < A < 
oo, independent from (C„) and from Z, and let 

J{t)=A,^. 

Remark 2. To define this process in terms of kernels, let -Pi"^i,T„_i he the 
exponential distribution with intensity A (regardless o/x„_i, t„_i Define Vt = 
max{n : t„ < t}. Besides, let tTxI^-i^t^, representing the increments at the jump 
instant Tn, be the normal distribution centered at {p — l)xn-i with variance 
(l-p')/2. 

We call such a process a Poisson-auto-regressive process (of order 1, covari- 
ance p and intensity A). 

Proposition 1. A Poisson-auto-regressive process J is wide-sense stationary. 

Besides, the random variable J{t) has a centered normal distribution with vari- 
ance 1/2 for every t. The covariance between J{t) and J{t + r) is given by 

gA(l-p)r/2. 

Proof. We compute the distribution and the covariances of J{t) conditioning 
on the number of jumps of J: 

oo 

P ( J(i) <x) = Y^ p,, (n)P ( A„ <x)=]P{Ao<x), 

n=0 

since An has centered normal distribution with variance 1/2 for each n, hence, 
so does J{t) for all t. In particular, J{t) is centered for all t. 
For the second moments we have: 

EJ{t)J{t-{-T) 

OO 

= Pu, (n)P {v{t, t + T\ = A;)E {J{t)J{t + r) | i/t = n, v{t, t + T\=k) 

n,k=0 

^ ^ f \j.\n ^ / \ \k k 

n,fe=0 n=0 ■ fe=0 
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where ^{t, t+r] denotes the number of jump instants in the interval (t, t +r] . □ 

Next theorem gives the mean number of continuous and discontinuous cross- 
ings through the level u by the process X on the interval [0, T]. For the sake of 
notational simplicity we consider only the up-crossings, but the case of down- 
crossings is completely analogous. 

Theorem 2. Let X = Z + with Z, J independent processes such that Z a 
stationary, centered Gaussian process with T{0) = 1/2 and paths and J a 
Poisson-auto-regressive process, then 

EUu^tJ^^{u) + \Tpp{u), 
V Ztt 

where ip stands for the standard Gaussian density function, A2 is the second 
spectral moment of Z, and pp(u) is the probability that a two dimensional, cen- 
tered, Gaussian vector with unit variances and covariance (1 + p)/2 belongs to 
the set (— oo,u) x (u, cx)). 

Proof. We begin considering the mean number of continuous crossings, as Z is 
Gaussian and stationary, we can apply Corollary [3l Besides, by Proposition [1] 
the density Px{t) — Px{o) = therefore 



Finally, an elementary computation shows that E Z{0) = •\/A7727r. This gives 
the first term of E Uu ■ 

Let us consider now the discontinuous up-crossings. The compensating mea- 
sure of the point process (t„, AA„)„, AAn — An — is Xdt F{dy), where F 
is the normal distribution centered at {p— 1)A^^_ and with variance (1 — /3^)/2, 
see [H eq. 4.64]. Hence 

EC/^ = E^ J l{X{t-) <u,X{t-) + y> u)\dtF{dy) 

= AE / dt I l{X{t~) <u,X{t~)+y>u}F{dy). 



Actually F is the distribution of A.A^^ conditioned on the random vector 
{Z{t),A,^_), so 

EU^ = XE [ P {X{t-) <u,X{t-) + AA^^ > u\ Z{t),A^^_) dt 
Jo 

= AE / V p,^ (n)P {Z{t) +A„<u, Z{t) + A„ + AA^+i > u \ Z{t), A„) dt 

•^0 n=0 

A / V P., (n)EP {Z{t) + A„ < u, Z{t) + An + AA„+i > u \ Z{t), A„) dt 

Jo n=0 

A / ^Pv, (n)P {Z{t) + An < U, Z{t) + An + AAn+l > u) dt 

Jo n 
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where we have conditioned on the number of jumps and used Fubini's theorem. 
By the stationarity of Z and J , the latter probability does not depend on t 

and n, hence 



E [/^ = AP (Z(0) + Ao < M, Z(0) + ^0 + A^li >u) \ V Vu, {n)dt 

n=0 

= ATP (Z(0) + Ao<u, Z{0) + Ao + AAi > u) . 

Now, observe that Aq + AAi = Ai = p^Q + \/l — p^^i and that the vector 
(Z(0) + Ao, Z{Q) + Ai) has centered normal distribution with variances one and 
covariance (1 + p)/2, therefore, the latter probability equals Pp{u). This gives 
the second term and completes the proof. □ 

Corollary 4. As u ^ oo, we have 

lim ^ = 0. 

Thus, the continuous crossings predominate, in mean, for high levels. 

Proof. Let (X, Y) be a two dimensional centered Gaussian vector with E X"^ = 
EF2 = 1/2 andEXF= (l + p)/2, then 

Pp{u) = p (X < w, r > w) < p (y > w). 

Then ^ 

Pp{u) < 1 - < -ip{u) = o((p(m)). 

The result follows. □ 

Remark 3. The processes Z and J contribute in 1/2 to the variance of X, 
more generally we can define 

X = aZ + Vl - a^J 

for a G [0, 1] and Z, J independent, centered, stationary with variance one. 
Hence, 



n 



= aT\l ^^2(^(0)) ^^^-^ = aTE \Z{Q)\ ^p{u) 



= aEN-{Z). 



Therefore, the mean number of continuous crossings of X is proportional to that 
of the continuous process Z , the constant of proportionality being a, that is, the 
standard deviation of the first summand aZ. 
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We now use Rice Formula to get upper and lower bounds for the tail of 
the distribution (i.e: for the overlife function) of the maximum of the process 
X with Gaussian stationary continuous part and Poisson-auto-regressive jump 
part. 

Thus, we consider the maximum random variable of X, 

M{T) = max{X(s) : < s < T}. 
These bounds are based on the elementary relation 

{M{T) >u}^ {X(0) > m} [+) {X(0) <u,Uu> 1}, 
where W denotes the disjoint union. It follows that 

P(M(T) > u) < P(X(0) >u)+¥{Uu> 1) < P(X(0) > u) +EC/„, (2) 

and 

P {M{T) >u)=¥ (X(0) > It) + P > 1) - P {Uu > 1, ^(0) > u) 

> P (X(0) > u) + E [/„ - C/„.[2] -FiUu> 1, X{0) > u), (3) 

where a^j = a(a — 1) is the Pochammer symbol and C/„.[2] = (CAt)[2] 

The following theorem contains the upper bound for the tail of the distribu- 
tion of the maximum of X on the interval [0, T]. 

Theorem 3. As u ^ oo the overlife probability of the maximum verifies 

P {M{T) >u)<l- $(m) +tJ ^ip{u) + \Tpp{u), (4) 

V ZTT 

where Pp{u) is defined in Theorem\^ 

Proof. Is a direct consequence of Theorem [2] and formula ([2]) . □ 

Remark 4. Furthermore, if we denote the r.h.s. of ^ byrhs{u), by Corollary 
[2] we have: 

rhs{u) ^ T \l ip{u). 
>oo V 27r 

The goal of the rest of this section is to show that this upper bound is sharp. 
In order to do that, we use the lower bound given in equation ^ for the overlife 
function of the maximum M{T). So we have to deal with the second moment 
of the number up-crossings. 

Theorem 4. Let the processes X and J be as in Corollary\3i If in addition 
J is a Poisson-auto-regressive process, and Z verifies that T{t) ^ ±1/2 for all 
T and the Geman condition: 

9'(t) 

— —dr converges at t — 0, (5) 
where is defined by T{t) = 1 — A2T^/2 + 0{t). Then 

P (M (0 > w) = 1 - $(u) +tJ^ Lp{u) + o{ip{u)) 

V 27r 
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Proof. It suffices to show that the additional terms in ([3]), w.r.t. ([2]), are o{(p{u)). 
Since 

Uu,m = u^iUu - 1) = c/^,[2] + c/,^,[2] + 2 ut 

the proof is divided in several steps, considering separately each one of the 
resulting terms. 

Claim 1. Wt have E,U^ — ^i'Pi'^))- 

Step 1. We begin with an upper bound for the second moment of C/„. 



</ / ^ E|Z^(s)Z-^(f) |X(s) = X(t) = u,i^,, =m,i^i_, =ji 

7n,n—0 

■ Px(s),x(t),^,.ut-s ™' n)dsdt. (6) 



I I ^\z^{s)Z+{t)\X{s) = X{t) 

Jo Jo „_n '- 



This bound is enough for our current purposes, but one can easily obtain the 
equality by the method of approximation by polygonals. 

We adapt the proof of the Rice formula for the factorial moments in [31 
Theorem 3.2]. Let C„ be the set of continuous up-crossings of X in [0,T], 
C2 = C„ X Cu and for any Borel set J in [0,T]2 let ^l{J) = 4t^{Cl n J). It 
follows that C^^ [2] — /^([0,r]^ \ A), where A is the diagonal, that is, A = 
{(s,t)e [0,T]2 '■.s = t}. 

Take Ji and J2 disjoint intervals in [0, T] and let J = Ji x J2(c [0, TJ^ \ A), 
then 

/l(J) = Uu{Jl) ■ Uu{J2) 

= ^"^.1^1 Z+{s)l{\X{s)-u\<5}ds- f Z+{t)l{\X{t)-u\<d}dt 
s^o {2Sy Jj^ Jj^ 

= lim J J Z+{s)Z+{t)l{\X{s) -u\<5, \X{t) - u\ < 5}dsdt, 

where we applied Kac Formula on each interval, and noted that for S small 
enough the quantity in the limit becomes constant, so we can use the same 
mute variable 6 in both limits. 

Now, we take expectation on both sides and apply Fatou's Lemma and 
Fubini's Theorem to pass the expectation inside the integral sign. Then, we 
condition on the number of jumps of the process in the intervals [0, s], [0, t] and 
on the values of the process X at these points, that is (if s < t): 

EZ+{s)Z+{t)l{\X{s) -u\<S, \X{t) -u\<S} = 

/ / X! ^ Z+{s)Z+{t) \ X{s) ^ x,X{t) ^ y,i^s ^ m,iyt ^ m + n 



m,Ti— 



Px(^s).x{t)\u,=-m,,^t^,^,=n{x, v)Pv, im)p„t [m + n)dxdy 



Observe that, under these conditions X{s) = Z{s) + A,n and X{t) — Z[t) + 
Am+n are jointly Gaussian, thus the conditional expectation is well defined 
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(and may be computed) via regression. Besides, this fact yields the regularity 
conditions needed for the integrand. In fact, the conditional expectation and 
the joint density function are continuous for t € [0, T] and x, j/ in a neighborhood 
of u. Hence, we can pass to the limit w.r.t. u inside the integral sign. 

So far we have stated that E fJ,{J) is bounded by the integral in the r.h.s. 
of ^ for any interval J C [0,T]^ \ A. As both sides represent measures on 
[0, T]^ \ A, the result follows by the standard arguments of Measure Theory. 

Step 2. Actually we can take inequality ^ a little further: 

^K.m < 2 /"^(T - r) E \z+z+ I ro(0) - r„(r) = u ■ 

■ PYo{o),Y„(T){u,u)Pur {n)dT, 

where we set Yk{t) = Z{t) + Ak, for t e [0, T] and fc e N. 

In effect, conditioned on Vs = i^t-s = n (if s < t; similarly on the other 
case) wehave Ar(s) = Z{s)+Am =: Ym{s) a.ndX{t) = Z{t)+Am+n —■ Ym+n{t)- 
It is easy to see that the vector (Ym{s), Ym+n{t)) is independent from ly and has 
centered normal distribution with variances 1 and covariance Tt-s + p'V2i in 
particular, this law does not depend on s,t but on the difference i — s, neither it 
depends on m. Therefore, the conditional expectation in inequality (|6]) reduces 
to that in the r.h.s. of the claimed bound. 

Then, we factorize the density function as: 



= PYo{o),y^{t){u, u)py, {m)p^^ (n). 

Clearly X^mPi^sC'^) — 1- Finally we make the change of variables (s,t) 
{s^T — t — s), and obtain the desired inequality. 

In the next two steps we bound each factor in the integrand. 

Step 3. A standard regression shows that if Vr = n, then: 

T'(t)u 

E(Z(0) I C) ^ -E{Z{t) I C) - ^ ' 



l-(r(r)+p'V2)' 
and 

var{Z{0) I C) = var(Z(T) | C) = A2 ^ ' 



l-(r(r)+p"/2)2 



where C = {A"o(0) = A'„(t) = u}. Therefore, using the inequalities a+5+ < 
(a + 6)2/4 and (a + bf < 2{a^ + b^) it follows that E(Z+(0)Z+(r) | C) < 
i(mr(Z(0) I C)+variZ{T) \ C)). Them 

E(^+(0)Z+(r) I C) < A2 ^'^^^ 



1 - (r(r) + p"/2)2- 



Note that for n > 1, as p < 1, there is no problem when r — ^ 0, in effect 
1 - (r(T) + p"/2) ^ 1/2 - p"/2 > 0. 
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Step 4. The vector (yb(0), y„(T)) is normally distributed with variances 1 
and covariance r(T) + /9"/2, therefore, the exponential in the density is 

exp 



1 + r(r) + p"/2 
For n > 1 we can bound p" < \p\ < 1, hence 



^^n^ 1 + r(.) + p"/2 1 - ^^n^ 1 + r(.) + H/2 1 = 

For ri > 1, replacing in the inequality ([B]) of Step 1 we have: 
EC/^.[2] <2T 



. A2(l-(r(r)+p"/2n-r-(ry ; ^ , , 

(l-(r(r)+pV2)2)3/2 "^P\ l + r(r) + H/2^''^ 



•^0 „=o 



<9T^ , .vp 



Vl-((l + H)/2)2 I (3+|p|)/2 



,2 



o((p(u)). 



The case n = 0, when there are no jumps in [0, t], is treated as in [3l Proposition 
4.2], in particular, we need Geman condition to ensure the convergence of the 
integral. 

Claim 2. We haveEU^ [3] ^ oifiu)). 

By the arguments in Corollary |4] it suffices to show that 

for some constant c and u large enough. We can write 

n=0 fc=0 

hence, making the product and taking expectation, we have: 

00 n — 1 
n=2 l=k<l 

■ P (Z(rfc) + Afc_i < w; Z{Tk) +Ak>u; Zin) + < u; ^(r^) + > u) 

00 n— 1 

<E E p.TWP(^(^fe) + ^fc-i <«;^(^^) + ^^>^i)- (7) 

n=2 fc<f=l 

Besides, 

P (Z(Tfc) + Ak-i < u- Z{n) + A,>u) 

ds / dtp^^^r,\uT=n{s,t)f{Z{s)+Ak-i<u-Z{t) + Ai>u) 

ds / dtpr,,re\uT=n{s,t)F{Z{0)+Ak-l<U;Z{t-s)+Ai>u), 
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in the latter equality we used the stationarity of the process Z. The vector 
(Z(0) Z(t— s)^ As) is centered Gaussian with variances 1/2 and covari- 

ance r(T) + /-'''+V2. It is easy to check that -r(T) - /-'=+V2 < 1/2(1 + 
which is the covariance of the vectors Z + S and Z+|p|S'+-\/l — p^V ■, therefore, 
by the Plackett-Slepian inequality, see [3] Section 2.1], we have 

P (Z(0)+Afc_i < w; Z(t-s)^At > u) < P (Z+S* < u; Z^p\S^^\- p^V > u). 
This bound does not depend on s,i, hence 

P (Z(Tfe) + Afe_i < u; Z(t£) ^^Ag>u) 

<PiZ + S <u:Z+\p\S+ ^/^~^V > u) dsj^ dtpr^^rAuT=n{s,t) 

= f{Z + S<u;Z+ \p\S + ^/l-p^V > u) 

= P|p|(u), 

since Tk,Ti \ vt — n- is concentrated on [0, T]^. Finally, replacing in the equation 
([7]) we have 

oo n-l (XT)'^ 
n=2 l=k<e 

and the result follows. 
Claim 3. We have EC/„ = 
In effect 

EC/„,[2] =E(C/^ + C/„^)(C/^ + C/„^-l) 

= EC/^^]+EC/„^[2]+2EC/^[/„^ 

< E C/^p] + E C^u,[2l + 2^Em^E{U^r 

= E ;^nj2] + E t/„',[2] + + ^ [2] + E C/„^), 

where we used Cauchy-Schwarz inequality. The first two terms in the r.h.s. are 
treated in the previous lemmas, under the square root sign the first factor is 
equivalent to (p{u) and the second one is o{ip{u)). Thus, EJ7„ [2] = o{ip{u)). 

Claim 4. We have F{X{0) > w, C/„ > 1) = 0{ip{{l + S)u)). 

We follow the proof of the analogue assertion in [Sj Proposition 4.2]. The 
key fact is that the distribution of the process remains unchanged under time 
reversal t^^T — t. 

In effect, let us condition on the number of jumps vt = n, then it is easy 
to check that the (conditional) distribution of (Ai, A2, . . . , yl„) | z^t = is 
the same as the distribution of (A„, . . . , Ai) \ vt — n. Besides, the 

distribution of ti , T2 , . . . , t„ | vt ~ n is that of an uniform (ordered) sample of 
size n, so it looks the same from and from T. Since the construction of the 
process J depends on these elements and there is no difference if we start at 
or at T the claim follows. 
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In conclusion, 




Taking into account (|1]), this completes the proof. 



□ 



3.2 Stationary Gaussian continuous process plus CPP 

In this example we consider a centered, stationary, Gaussian process Z with 
r(0) = 1 and paths and an independent Compound Poisson Process (CPP), 
J with finite intensity A and standard Gaussian jumps. 

The process J is defined by the kernels -Pi"^i,T„_i, the exponential distribu- 
tion of intensity A > and 7ri"-*_j^_T„, representing the increments, the standard 
normal distribution. We can also define J in terms of random variables as fol- 
low, let (Tn)n&i be independent random variables with exponential distribution 
with intensity < A and {t,n)n&i be independent standard Gaussian random 
variables independent from (t„); then for f > define: 



Thus, J(i) is the sum of a random number of standard normal random variables. 

It is well known, see [9] for instance, that the compensating measure of the 
CPP is \dt^{dx), where $ is the standard Gaussian distribution. In particular, 
the compensating measure is deterministic. 

Denote by the centered Gaussian density with variance n, in particular, 

= </? is the standard normal density. Note that (p„ * = </Jn+m, where * 
stands for the convolution. 

Next theorem gives the mean number of up-crossings through u hy X. The 
mean number of down-crossings is analogue. 

Theorem 5. For the process defined above we have 



where X2 is the second spectral moment of Z and p — '^'^^iPn^n with p., 



Remark 5. In spite of the notation, the first term, which corresponds to the 
continuous up-crossings, does depend on A through the density function p, which 
has expectation and variance XT/2. 

Remark 6. Observe that the mean number of continuous crossings, the first 
term, has the usual form, see Remark]^ 





^^{vT > n}. 
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Proof. We begin with the continuous crossings. Using the formula in Remark 
[U we have 



EU^^ dt E [Z+{t) I Zit) = v]pzit){v)p.m)iu-v)d 



= EZ+(0) / dt Pz(o){v)pj(t){u~ v)dv, 

Jo J -oc 

where we used the stationarity of Z as in Coroharyll] RecaU that EZ+(0) 
■y/ X2/2tt. Now, by Lemma [U we decompose pj{t), so 



Ei 



n=0 

2^ Jo 



^u^Xtt- dty^P^tin) pz(o){v)Vn{u-v)dv 

V -/O . n J-oo 

oo 

(0) * V^n )(") 

n=0 

This integral is computed in the second item of Lemma [51 we get 
E C^u = r ^ — P {;.T > n + 1} ^„+i(u) 



fx /A 



n=0 



n=l 



= T a/ ^ 

V ZTT 



which gives the first term of the result. 

Now, we turn to the discontinuous crossings 



EUf^=E / l{X{t-) <u;X{r)+y>u}L{dt,dy) 

Jo J-oo 

= e/ / l{X{t-) <u:X{t-) + y > u}Xdt<S>{dy), 

Jo J -oo 

where, as we said before, $ is the distribution of the jumps. Since the compen- 
sating measure is deterministic we have 

pT poo 

^U^ = >^ / E1{X {t-) <u;X{t-) + y>u}df^>{dy) 

Jo J -oo 

{X{t-) < u;X{t-) + y> u)dt^{dy) 



J -oo 



= \ F{X{t-) <u;X{t-)+^> u)dt, 
Jo 
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being ^ a standard normal variable independent from X{t~). Now, we condition 

on X{t~), since, for fixed, t almost surely J{t~) = J{t) we may use the same 
computations as in the computation of E A^°: 

EU^^xf [ V{^>u-x\X{t-) = x)px{t-)(x)dxdt 

Jo J-oo 

/u _ i-T 
^{u — x)dx / Px{t-){x)dt 
-oo Jo 

= XT f $(u - x)p{x)dx. 

J — oo 

This completes the proof. □ 

Next, we compare the mean numbers of continuous and discontinuous up- 
crossings through the level u by as u ^ oo. 

Corollary 5. As u ^ oo the mean numbers of continuous and discontinuous 
up-crossings through the level u are of the same order. More precisely 

X /Ott 

lim -—Em< lim EU^< lim X\ —Em 

u^oo a/ A2 It— >oo >oo y ^2 

Proof. We bound from above and from below Et/^. First, we have 

/u poo 
- y)piy)dy = / ^{y)p{u - y)dy 
-00 Jo 

p2u POO 

= / ^{y)p{u - y)dy + / <^{y)p{u - y)dy 

Jo J2u 

p2u f-2u 

> / '^{y)p{u - y)dy > p{u) / '^{y)dy, 
Jo Jo 

where we used that p is even and decreasing on [0, 00). Furthermore 

/•2u POO 
lim / ^y)dy= / ^y)dy = E^+ = J- 

f°° 

EUi>XTp{u) ^{y)dy ^ ^p{u). 

Jo v27r 



Therefore 



On the other hand, 



00 />00 ^ PQC 

EU^ = XTj^Pn / Hu-y)'Pniy)dy < xrJ^Pn / ^ 

n=l n=l -^-"o 



{u-y)ipn{y)dy. 



Note that this is a convolution formula, for the overlife function, of two (inde- 
pendent) random variables, say Z ~ $ and V„ ~ (p„. Then we can write the 
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latter integral as P(Z + l/„ > u). Furthermore, since $ is the Gaussian stan- 
dard distribution and Lpn is the Gaussian density with zero mean and variance 
71, this probability equals V {Vn+i > u) — <^{u/\/n + 1). Thus 



where 5{u) := J2^=iPni''^ + l)</'n+i(ti) — p(u). In order to obtain the desired 
result, it suffices to prove that S{u) a.s u 00. In effect 



\S{u)\< 



[Pnin + 1) -Pn+l] Vn+liu) 



n=l 



+ \Pi'Piiu)\ 



It is clear that the second term of the r.h.s tends to zero when u — > oo. Let 
us look at the first term, note that J2'^=iPn(j>' + 1) is, roughly speaking, the 
second factorial moment of the Poisson distribution, in effect 



^ p„(n + 1) = 1 + — ^ nP (i/T > = 1 + — ^ n ^ p.^ (fc) 

n=l n—1 n—lk—n 

^ 00 A; ^00 



k=ln=l k=l 

Thus, the first term is a mixture of Gaussian densities (times a constant), 
hence it tends to zero. Then S{u) — ^ as u — > 00 and 

EU^<XT{p{u) + S{u)) - XTp{u). 

Now, the result follows joining the two obtained bounds for E and using 
Theorem El □ 

Remark 7. Note that in this case the discontinuous crossings through the level 
u are not negligible w.r.t. the continuous crossings when u tends to infinity, in 
contrast with the situation in the example of the previous section. 

With the usual procedure, i.e. equation ©, we obtain the following result. 

Corollary 6. As u — > 00 the tail of the distribution of the maximum is bounded 
from above by 

P {M{T) >u)<V (X(0) >u) + T \l^p{u) + \j Hu~ y)p{y)dy. 

4 Generalization of Borovkov-Last's formula 

Borovkov and Last in are interested in the continuous crossings through a 
level u by a stationary Piecewise Deterministic Markov Process. A process X 
of this kind, starts at a random position, then jumps a random quantity at 
random times but moves deterministically between jumps. 
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Such a process is described by a general point process (t„, ^„)„, as described 
in Section [5J and a (non-random) rate function /i : M — R. More precisely, the 
process X has jumps at the points {r„), the magnitude of the junp is ^„ and on 
the interval [t„, t„+i), X{t) follows the integral curve of /i with initial condition 
X(r„)=X(r„-)+e„. 

Note that the jump part of the process is not independent of the continuous 
one. 

Let Dfj, = {u : fi{u) = 0}. Observe that if //(u) > (resp. <), the continuous 
crossings of level u can only be up-crossings (resp. down-crossings) . 

Next theorem extends Borovkov-Last Formula to the non-stationary case. 

Theorem 6. Let u ^ and assume that fi and px(t) ^'^e continuous w.r.t. x 
in a neighborhood of u and t £ [0,r]. Then: 

E7V^ = |/i(u)| / px(t){u)dt 
Jo 

Proof. For the levels u ^ we can apply Kac counting formula pathwise for 
almost all paths of X. In effect, the continuity of fi implies that the paths are 
of class between the jumps and that X{t) — fi{X{t)) for almost all t G [0, T]. 
Since X{t) has a density, the value u is not taken at the extremes of the interval 
neither at the jump points almost surely. Furthermore, by the continuity of this 
density, there are not tangencies at level u. 

Now we take expectation on both sides of Kac Counting Formula and ob- 
serve that the number of continuous crossings of the level u ^ is bounded by 
the number of jumps -1-1 of in [0, T] . In effect, the sign of /i(u) determines the 
direction of the continuous crossings of m, so, between two continuous crossings 
must be a discontinuous one in the opposite direction, thus there is at most one 
continuous crossing at each one of the intervals of the partition tq, . . . , r^^ , T. 
Then, since vt in integrable, we may pass to the limit under the expectation 
sign: 



1 '■^ 



EiV,'; ^ lim— / E 
" Sin 25 Jo 



\^it)\'i-{\X(t)-u\<S} 



dt. 



Now, X{t) is a deterministic function of X{t), namely X(t) = iJ.(X(t)), so the 
integrand is simply the expectation of a function of X{t), therefore 

EK = l|m^^ H^)\Px{t)ix)dxdt. 

By the continuity of the integrand and the compactness of the domain we can 
pass the limit inside the integral w.r.t. t. Then, the result follows by the mean 
value theorem. □ 

As a corollary, when X is stationary, we obtain Borovkov-Last's Formula. 

Corollary 7. // in addition to the conditions of Theorem the process X is 
stationary, then 

where p = px(o)- 
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Remark 8. Note that on this case the (net) number of discontinuous crossings, 
namely the difference of down and up crossings D'^ — U^, is related to the number 
of continuous crossings, actually we have: 

EDt-EU^ = EN^ + sgni^I{u)){FiXiT)<u)-P{X{0)<u)) 

5 Proofs 

First, we prepare some auxiliary results. 

Lemma 1. Let T, X,Y, Z be random variables, then 

/oc 
E{T\X^x,Y^y,Z^ z)pz\x=x,Y=v{z)dz. 
-OO 

and 

E{T\X^x,Z^z)px\z=z{x) 

E{T \ X ^ x,Y = y,Z = z)px.Y\z=z{x, y)dy. 



Proof. These equalities follow directly from the properties of conditional expec- 
tation, see [121 Page 215]. For example, for the first one take F2 = {X,Y) and 
F^ = {X,Y,Z). □ 

Lemma 2 (Auxiliaries). Consider a CPP {J{t)) with standard Gaussian jumps, 
then 

1. the density of the CPP can be written as: 

00 

Pj{t){x) ^p^t{n)ipn{x). 

2. The integral 

pMdt = [UT > n + 1) - ip (r„+i < T) 

3- IfPn - wF(!yT > n), then Y.n=iPn = 1 
As usual we denote the time epochs by r„ and the number of jumps in [0,i] by 

Proof. 1. Conditioning on the value of vt we have: 

00 

Fjit) (x) = P {J{t) <x) = Y, P^t {nW (Jit) <x\i^t^n). 

71=0 

Now, if i^t = n, J{t) is the sum of n independent standard Gaussian random 
variables, hence, the conditional probability in the r.h.s. of the latter equation 
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is the distribution of centered normal random variable with variance n. The 
result follows taking derivatives on both sides. 

2. By definition, p,^^ (n), as a function of t, is equal to the density function of 
the Gamma distribution, with parameters A, n + I, divided by A, furthermore, 
it is well known, |15) , that this is the distribution of t^+i, hence 

PvAn)dt = . 

The result follows since the events {t„+i < T} and {vt > n + 1} coincide. 

3. It follows directly from the facts that vt has Poisson distribution with 
mean AT and that for a non- negative integer valued random variable X: ¥,X ^ 

Proof of Theorem [TJ Let us start with the formula for the mean number of 
continuous crossings, we compute the expectation by conditioning on the paths 
of the pure-jump part and use the proof of the non-Gaussian case on [^. 

Then, we condition on the number of jumps, vt — n, on the jump instants, 
Tfe = tk, thus 

EN-^=E (E [N^ (A-, [0, T]) \i^T^n;T = t]) 

OO „ 

= y^p^An) / Px(t) E [N^ (A-, [0,r]) I = n;T = t] dt, 

where we set t = (ri, . . . , t„) and t = . . . , 

Now, we look at the integrand and, since the number of crossings is additive 
w.r.t. the interval, we split the interval [0,r] as the union of the intervals 
Ife [tk-i,tk), then 

n 

E [N^ {X, [0,r]) I = n;T = t] = ^E [7V^ {X,Ik) \ i^t = n; r = t] 

fc=i 

Each term can be written, see Lemma ([1]), as 




IE [N^ {X,Ik) \vT ^n;T ^ t;J{tk-i) = y]pj{t^,_i)\i,T=n-i-=t{y)dy. 



Now, conditionally on vt = n;T — t and J{Tk-i) = y, the process X can be 
written as -Z -|- y on . Since Z verifies the conditions Al, A2, A3 and AA on 
Xk so does the process 2 + y, therefore we may apply Rice Formula, see [3], on 
each interval under these conditions to obtain 



E [N^ {X,Ik) \yT=n-T = t; J{tk-i) = y] 



E 



\Z{t)\ I X{t) ^u,UT^n-T = t, J{Tk-i) - y 

■ PX{t)\uT=n;T=t,J{Tk^i)=y{u)dt 

Each one of these expressions should be replaced in the previous one. 
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Finally, we have to integrate (the conditions) , for notational simplicity, let us 



write g{n,t,y) 



\Z{t)\ I X{t) ^u,VT^n-T^t, J(Tfe_i) = y 



E \Z(t)\\ X{t) u,VT = n:T = t and = E \Z{t)\ \ X{t) = u,vt ^ n 

Let us perform the integrals one by one, starting w.r.t. y, (use Fubini) then 

9[n,t,y)px(t)\uT=n-,'r=t,J{Tk)=y{u)pj{tk)\,,T=n-,-r=t[y)dy 

g{n, t, y)p(x{t),j(r^,))\vT=n--r=t{u, y)dy 



'l{n,t)px{t)\ 



UT—Tl\T—t 



Where we used Lemma ([T]). Now, we sum the integrals over k and integrate 
w.r.t. t: 



dt 



[0,T]" 



9{n,t)px(t)\vT=n--r=t{u)p-r\vT=n{^)dt 



dt 



[0,T]" 



gin-, t)p^x(t),-r)\vT=n{u-t)dt 



9{n)px(t)\uT=n{u)dt. 



By the same arguments one can remove the condition on vt- The result follows. 

Now we proceed to the formula for the mean number of discontinuous cross- 
ings through level u. For the definitions used below see [9]. 

Clearly, X only can have a discontinuous crossing through u at the points 
T„; n — \, . . . ^ Vt, the jump of X at each one of these points is due to the jump 
of J'. Hence, we consider the Marked Point Process {{Tk,^k) : k > 0) associated 
to J' on [0,00) X M, which defines a Random Counting Measure fi{dt,dy), in 
terms of which we can write: 



k=l 



0<t<T 



Nt = J2 mX{r^) u){X{r^) + AX(r,.) - u) < 0} 

= ^ l{{X{t-)-u){X{t-)+U^u)<Q} 

l{{X{t-) ~ u){X{t-) + y~u)< 0}M(dt, dy) 

[0,T]xR 

Is easy to see that this RCM has a compensating measure denoted by L{dt, dy), 
see [9] again. Taking expectations on both sides we have: 



E< = E [E(7V„^|Z)] 



E 



E 



[0,T]xR 

= E 



l{{X{t-) - u){X{t-) +y-u)< 0}f,{dt, dy) \ Z 

l{{X{t-) - u){X{t-) +y~u) < 0}L{dt, dy) 

[0,T] xR 
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where we used that, conditioned on the integral is done w.r.t. the RCM [i 
associated with J which coincides with the integral w.r.t. the compensating 
measure L{dt,dy) since the integrand is predictable, in fact it is a function of 
t~ . Finally we integrate with respect to Z and the result follows. □ 
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